UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRARY ROUTINE V- 5 - 18^ 



TITLE 

TYPE 

NUMBER OF WORDS 

TEMPORARY STORAGE 
PRESET PARAMETER 



DURATION 



ACCURACY 



Spherical Bess el Functions 

Closed 

59 (38 plus the 21 of Sine-Cosine Routine T - 5 

■which is attached . ) 

0-6 inclusive and ( k + 1 ) locations at S3 

e 

S3 contains the address at which the table of 

functions j (x), ^(x), ».. , <j' k (x) is to 

be placed* 

(l8 + 2*5 k ) as where k is the order of the 

highest Bessel Function to be found. 

Extremely variable and dictates that this program 

shall be used with care* See discussion under 

"Method and Error Discussion" below. 



DESCRIPTION AND ENTRY 



-n 



Enter with the scaled argument 2 x in the 
accumulator (A) and the program parameter: 



p 


50 (l6k + n)F 
50 pF 


p + 1 


26 qF 



1 < k < 15 

— e — 



NOTE 



**n 
where 2 is the scale factor on the argument, 

k is the maximum order of the Bessel Functions 
e 

to be found, and q is the address of this program. 

On exit j (x) will be found in rS3 i*e« j (x) in S3* 

j (x) in 1S3, J 2 (x) in 2S3, . . . , j k (x) in k e S3« 

e 
This program contains Sine-Cosine Library Routine 

T 5 - 157 as an appended part. It is available for 

use in other connections, so it need not be placed 

a second time in the machine. Merely enter it at 

word 38L of this program- 
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METHOD AND ERROE DISCUSSION 



The Spherical Bessel Functions j-^fx) are defined 
in terms of the ordinary Bessel Functions 



J k (x) ^ 



J k* x) =^£F— J k + l/ 2 W ' ( x ) 



Explicitly : 



sin x 



Jo (x) = x ' (2) 



j ± (x) - _£i|_i . _£2iL2£ , ( 3) 

X X 

j 2 (x) = / — ^—* — ) sin x ^- cos x, 

V x x/ x 



etc. 
They satisfy the recurrence relation; 



W x) = {2k + x) V x) - J k -i (x) - w 



Although the argument of j, (x) need not he less 

than 1 and hence in scale, |j (x) | < 1 and 

k 
hence is in scale. 



In this routine the sine-cousine routine 
is used to evaluate J n fcO and j,(x) using (2) and 
(3)» The relation (k) is then used to 
construct the table of functions . 
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Although this method is simple from 
a programming point of view, since it 
obviates the tedious evaluation of series, 
it can "be subject to large numerical error. 
This error in j (x) may be as large as of 
the order of: 

(2k ~ l)t! x~ k 2~ 57 (5) 

when x < 2k + 1, where 
(2k - 1)51 Z (2k + 1) (2k ~ 3)-.- (5) (3) (1) 

When x > 2k + 1,- the error is of the order of 

-38 
k2 * These estimates give the order of 

magnitude of the upper limits of the error. 

Table 1 gives the error as determined 

experimentally from consulting tables of 

these functions . This table could not be 

11 
read to better than 1 part in 10 

For x > 10 and j < 13, no error could 
~ k ~ 
be detected using the number of significant 

figures available in existing tables « 

All of these tests were made with an 

optimum choice of the scaling factor 

2~ n s i.e. n as small as possible consistent 

-n . 
with 2 x remaining m scale. 

For a given argument , in the regions 
of k labeled "worthless" the answers lose 
significance to such an extent that the pro- 
gram starts to be subject to division hangup 
in the later loops in attempting to compute 
the higher j • 

It should be noted that in the regions 
labeled "worthless" for a given argument, a 
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good approximation to j is zero* This 
approximation will (in these regions ) be in 
error by less than the entry just above 
the region. This is because (in order of 
magnitude): 

.8 



■k 



'k 



'k 



(.01) 


< 


10 


for k > 3 


(-06) 


< 


io- 8 


for k > k 


(0.1) 


< 


10" 7 


for k > k 


(1.0) 


< 


io-T 


for k > 8 


(1.5) 


< 


io" 7 


for k > 9 


(3*0) 


< 


10" T 


for k > 12 



DAT E September Ik, 1933 

2. HujrvHMsfitf 



PROGRAMMED BY' 
APPROVED BY 







dh * 



mge 



TABLE 1 



Error in j. For Argument 



Function 



.01 



10 



-10 



10 



■""■"■p 



*o6 



10 



-11 



5 x 10 



-10 



10 



2 x 10 



worthless! 2 x 10 



36" 



0.1 



10 



-11 



1.0 



2 x 10 



-10 



5 x io 



-9 



3 x 10 



< 10 



-11 



< 10 



-11 



< 10 



-11 



10 



-11 



1-5 



< 10 



-11 



< 10 



-11 



< 10 



-11 



10 



-11 



5*0 



< 10 



-11 



< 10 



-11 



< 10 



-11 



< 10 



-11 



5.0 



< 10 



-11 



< 10 



-11 



< 10 



-11 



< 10 



-11 



10.0 



< 10 



-11 



< 10 



-11 



< 10 



-11 



< 10 



-11 



h 



worthless 



worthless 



7 x 10 



-11 



5 x 10 



-11 



< 10 



-11 



< 10 






< 10 



-11 



6 x 10 



-10 



k x 10 



-10 



< 10 



-11 



< 10 



-11 



< 10 



10 



2 x 10 



-9 



10 



-li 



10 



-11 



< 10 



-11 
■11 



6 x 10" 



3 x 10 



-9 



2 x 10 



-11 



2 x 10 



-11 



< 10 



-10 



worthies g i 8 x 10 



8 x 10 



-11 



7 x 10" 13 < 10 



-10 



i 

VJl 

i 



-1 -1 (mm 

worthless k x 10 2 x 10 1 < 10 



•10 



'10 



'11 



2 x 10 



-9 



8 x 10 i < 10 



-9 



2 x 10 



3 x 10~ 9 ■< 10 



-9 



'12 



ii 



Y 



V 



worthless 



V 



v 



v 



-8 



10 



6 x 10' 



< 10 



< 10 



< 
r 

H 

CD 



LOCATION 



4 



10 



11 



ORDER 



40 3F 
K5 F 
42 25L 
10 24f 
42 10L 
01 4F 
k2 6L 
42 % 
42 16L 
42 19L 
42 28L 
50 5F 
75 37L 
00 (n)F 
4o 4f 
50 7L 
26 38L 
00 IF 
40 5F 
10 (n)F 

66 3F 
S5 (k e )F 

50 5F 
40 S3 



NOTES 
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By 3 



By 3 



By 2 





Plant Link 



Plant k 



n In R- (EHA) 



Plant n 



2~\ in R, 



R, 



2 ~* (JL.) m H. 

x it ■ 1 

( X / ) mod 2 In E. 



X 



(-JL.) mo d 2 -» 4 



2" sin x in R. 



sin x In R- 



sin x -» 5 

*-n 
2 sin x in R_ R 

1 2 



V5-184 



J o 


in R 2 








J o 


in R* 


(End 


const, 


) 


-n . _ 
x m R p 








3 


-*S3 









V5-18^ 



LOCATION 



12 



13 



Ik 



15 



16 



IT 



18 



19 



20 



21 



22 



23 



ORDER 



TJ 3F 


Uo 6f 


LJ 


4f 


50 


13L 


26 38L 


00 


IF 


40 ^F 


ks 


31L 


L5 


5F 


10 


(n) F 


^4-0 




S5 


F 


50 


3F 


TO 


4f 


lA 


F 


10 


(n)F 


66 6f 


S5 


F 


^0 


1S3 


L5 


11L 


k2 


32L 


f4 56L 


42 


29L 


f4 36L 



NOTES 
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By 4 



By k 



From 35 



x in R n 

-2n . c 
x in o 



— + l/2 in R n 

it ' 1 



2" 1 


cos 


x in 


R i 


Cos 


X 


in R 




Cos 


X 


— » 


k 




Clear 


k 


counter 


s^-n 


X 


inR 1 




2 ~ 


sin 


X 





-n 



(2~ Sin x) (insp) -^ 



~n 



(2 Sin x) (ksp) -> R. 



- 2 x cos x in R_ R„ 

2*" (sin x * x cos x) is L R~ 

2" (sin x - x cos x) in B, R £ 
j in R 2 

J l in R l 

address S3, 1&3, and 2S3 




LOCATION 


ORDER 


2k 


42 3^L 




F5 31L 


25 


LO 10L 




32 ( )F 


26 


L^ 10L 




ko 31L 


27 


50 36L 




00 IF 


28 


Fk 36L 




10 (11 )F 


29 


10 4f 




75 (1S3)F 


30 


00 3F 




36 31L 


31 


66 3F 




S5 ( )F 


32 


hO F 




LI (S3)F 


33 


10 IF 




Lk F 


3^ 


00 IF 




kO (2S3)F 


35 


F5 32L 




26 22L 


36 


00 F 




00 F 



By 1 



By 5 



By 23 



By 22 



By 2k 



J 



NOTES 
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IV5-181* 



-39 

k 2 ^ in R 

(k ~ k ) 2" 39 m R 

© Done - out by link 

0Not done 

k x 2 ~ 59 in R IT 

Clear R 

2k 2*' 59 

(2k + 1) 2*" 59 

(2k + 1) 2" n in R 2 

(2k + 1) 2" n " 4 in R 2 

(2k + 1) j 1 (2" n ' 4 ) in R x R 2 

(2k + 1) j x 2~ n ~ 1 in R x R 2 



Waste 



(2k + 1) 



2x 



in R, 



and in R.. (Counter) 
(2k + 1) 



a 



2x 



in 



- J 







-1 



-h 2 

•2/ 
J / 2 

j 2 in R x 

J 2 -» 2S 3 

Step addresses "by 1 

and repeat 

- 



LOCATION 


ORDER 


NOTES PAGE k 


37 


00 F 
00 318 309 886 l84j 




- 1 


38 


50 39L 


- 


~) 




26 999F 




j Interlude to 


39 


00 F 
00 38L 




r~ 


ko 


00 F 

26 38L 






kl 


26 IN 


1 




k2 


Routine T5 (Sine Routine) 




ends this tape. It will 




"be loaded at 38L* 



V5-l8ij- 



JNS, 
DH ' 



:mge 



